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I. INTRODUCTION 

Quantum control of many-body states is essential for 
new quantum technology applications. Ultracold atoms 
trapped in optical lattices provide an excellent system 
to probe the coherent dynamical control of multiparticle 
states. A paradigm of the coherent control of cold atoms 
in optical lattices was provided by the superfluid (SF) 
to Mott-insulator (MI) transition via adiabatic ramping 
of the lattice potential More recently, a theoretical 
proposal |2|] to drive a SF-MI transition by using a fast 
oscillatory force has been experimentally demonstrated 
0. In this work we analyze the role played by disorder 
for interacting atoms in a lattice in the presence of a fast 
non-resonant oscillatory potential. 

The coherent dynamics of tight-binding systems in the 
presence of external fields is a very timely subject and 
includes phenomena such as Bloch oscillations, Wannier- 
Stark localization, caused by static fields, or dynamical 
localization due to the interaction with ac fields 0, 
The energy eigenvectors of a quantum particle in a peri- 
odic lattice are Bloch waves. The corresponding eigenval- 
ues form energy bands, and a single particle is delocalized 
over the lattice. Dunlap and Kenkre [6] showed that an 
external ac field localizes the particle when the ratio of 
the field magnitude and the field frequency is a root of 
the ordinary Bessel function of order zero. The use of 
Floquet theory gives rise to a great conceptual simpli- 
fication [1]; the structure of the quasienergy bands de- 
pends on the parameters of the time-dependent external 
perturbation and therefore by fine tuning the amplitude 
and frequency of the field, tunneling and the consequent 
localization phenomena can be manipulated. The dy- 
namical suppression of tunneling persists in the presence 
of interactions Q and it was proposed to induce localiza- 
tion for a single particle in the presence of disorder f^. 
Here, we combine the fast oscillatory force, the interac- 
tions between atoms, and disorder for a one-dimensional 
tight-binding system. 

The interplay between disorder-induced localization 
and interactions is one of the main open issues in con- 
densed matter physics. In the context of bosonic atoms 
trapped in optical lattices it has risen much interest in 
recent years 0, Q and the exact ground state phase dia- 
gram of the disordered Bose-Hubbard (BH) model is still 



an open question [TO, [O, [lH, [TO, [13] • For commensurate 
fillings and strong interactions in the absence of disorder, 
there exist insulating MI states characterized by an inte- 
ger fixed number of bosons per lattice site and a gap in 
the excitation spectrum. In the presence of disorder an 
additional insulating phase appears, the Bose-Glass (BG) 
phase characterized by a gapless excitation spectrum and 
finite compressibility. For weak interactions and strong 
disorder an Anderson-Glass (AG) may appear where one 
or more condensates localize in Anderson states. 

It has been shown previously that adiabatic ramping of 
a fast oscillatory force can drive an ordered system with 
commensurate filling from a superfiuid regime into the 
localized regime as it effectively reduces and even sup- 
presses the tunneling between neighboring sites [H . This 
force periodic in time but constant in space can be pro- 
duced in the laboratory through an effective acceleration 
induced by a linear-in-time dephasing of the lasers form- 
ing the optical lattice [s', . By changing appropriately 
the parameters of this acceleration, its amplitude and 
frequency, it becomes possible to control the transitions 
between different states of the BH Hamiltonian. We show 
here that adiabatic quantum evolution towards a local- 
ized state is still possible in the presence of disorder both 
for commensurate and non-commensurate fillings. The 
structure of the quasienergy levels becomes more com- 
plex when disorder is introduced. The breakdown of de- 
generacies eliminates spectral gaps giving rise to a denser 
spectrum, and in principle adiabatic evolution could be 
spoiled at least for reasonable time scales. 

In section|II]we introduce the Hamiltonian and the Flo- 
quet theory used to describe the dynamics of the system. 
Section IIIII discusses the time scales needed for inducing 
a transition into localized states for different ratios of the 
interaction and disorder strength. For commensurate fill- 
ing we find that in spite of the breaking of symmetries 
induced by disorder, the time scale needed for adiabatic 
evolution is similar for an ordered and a disordered sys- 
tem. Localization is not possible for an ordered system 
with non-commensurate filling via the adiabatic ramp- 
ing of the fast oscillatory potential [16] . We show that 
the breaking of symmetries induced by the disorder po- 
tential makes localization possible even at similar time 
scales than for the commensurate system. We also show 
how the adiabatic time scale shortens for increasing dis- 
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order strengths. Finally, we discuss the implications of 
our exact small system calculations for systems of exper- 
imentally relevant sizes. In Section [TVl we summarize our 
main conclusions. 



II. TIME-DEPENDENT DISORDERED 
BOSE-HUBBARD MODEL 

We consider the dynamics of N atoms trapped in the 
lowest Bloch band of an optical lattice at sufficiently low 
temperatures such that it can be described by the BH 
model. The corresponding Hamiltonian is (taking h = 1 
throughout) Q 
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where aj is the bosonic destruction operator for an atom 
localized in lattice site j = — M/2, .., M/2 and rij = a^^aj. 
The parameter J is the tunneling matrix element, U 
is the on-site repulsive interaction, while ej G [—1,1] 
denotes a site-dependent random external potential of 
strength W and C accounts for the external parabolic 
trap. Depending on the ratios W/U and U/ J the ground 
states of Ho show very different properties and, for large 
enough A^, different parameter values span the phase- 
diagram of the multiparticle system. 

For W = 0^ one can span the SF-MI quantum phase 
transition by changing the ratio U/ J. In the limit 
U/J-^0 the tunneling dominates and the system is SF, 
while for [// J ^ 1 the repulsive interactions dominate 
and the system is an insulator, with a gapped spectrum. 
For W and U/J^l the system is an insulator, 
but depending on the ratio W/U and the mean occu- 
pation number N/L (L = M + 1 is the total length of 
the system, always odd in our calculations), the ground 
state is either a BG or a MI. Both states are localized, 
in the sense that there is no transport through the lat- 
tice, and thus the density matrix has exponentially decay- 
in g off -diagonal terms. What distinguishes both phases 

[l2[ is the structure of the spectrum, which is gapped 
for a MI while it is not for a BG. When disorder dom- 
inates W/U > 1 the ground state is a BG, while for 
W/U < 1 different theoretical methods predict a MI or 
a BG ground insulator state depending on the mean oc- 
cupation number [TO, [O, [lH, [13] . For small disorder, in 
the regime U/J 0, the system is SF with long-range 
coherence while in the region U/W <C 1, the system con- 
denses to an Anderson state [I3 localized in the regions 
of strong disorder. These properties do not smear out in 
the presence of a shallow trapping potential. 

We impose a linear potential V{t) across the lattice 



which oscillates in time [sl, [l5[ and the total Hamiltonian 

H{t) = Ho^V{t), (2) 

V{t) = V cos(cjt + (5) ^ juj (3) 
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is then periodic with period T = 27r/a;. The Floquet 
theorem establishes that, due to the periodicity of the 
Hamiltonian, the solutions to the Schrodinger equation 
I ^(t)) = H(t) I ^(t)) are linear combinations of func- 
tions of the form 

\^%t))=e-'''\r{t)), (4) 

where \(t)''{t)) = |(/)^(t + T)). Inserting Eq. g]) into the 
Schrodinger equation one obtains the eigenvalue equation 

H^{t)\<p%t))=e\<p%t)), (5) 

where the Floquet Hamiltonian is defined as 

H^^H{t)-i^^, (6) 

and e is thequasienergy of the system. As pointed out 
by Sambe [l8|, Eq. ([5|) can be solved using the standard 
techniques developed for time-independent Hamiltoni- 
ans, provided we extend the Hilbert space to include the 
space of time-periodic functions. A suitable basis for this 
extended Hilbert space 7^ = H T is {| a) 1 n)}, where 
{\a)} is a basis for the Hilbert space H of the system, 
and we define {t\n) = e*^^^ with n integer, as the basis 
functions in the vector space T of time periodic func- 
tions. In this basis, Eq. ([5|) becomes a matrix eigenvalue 
equation of infinite dimension with an infinite number of 
eigenvalues. It can be proven that if is an eigenvalue 
with corresponding eigenvector |0^"(t)), then -\- poo 
with p integer is also an eigenvalue with corresponding 
eigenvector | 0^«+^^(t))= e^^^^ | (t)). Accordingly, the 
eigenstates \ (j)^'^'^^^ {t)) and \(j)^'^{t)) describe the same 
state in H apart from a global phase. Thus, in order 
to describe the time evolution in the Hilbert space we 
need only to solve the Floquet eigenvalue equation for 
£^ = Ca G (— ^cj, called first Brillouin zone in anal- 
ogy to the theory of periodic crystals. It can be shown 
that there are D distinct quasienergy states in the first 
Brillouin zone, if the Hilbert space H is D-dimensional. 
The set of cyclic wave packets | ^^°^(t)) in Eq. (|4]) ob- 
tained from the Floquet eigenstates provides a complete 
orthogonal basis in 7i at each t. The solution of the 
Schrodinger equation is then 

\^{t)) = U{t,to;V)\^{to)), (7) 

where 

U{t,to;V) = J2\<i>'-{t)){^'-{to)\. (8) 

a 

From now on the term cyclic state refers to |^^"(to)) 
with to being the beginning of a period. 
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We are interested in a regime in which uj is smaher 
than the energy separation to higher energy bands, so 
that our single-band BH model of Eq.([Tj) is still valid, 
but it is larger than the rest of the relevant energy scales 
of our system, i.e. a; ^ [/, J, W . In this high-frequency 
regime, one can show than the effect of the time-periodic 
external potential in Eq. ([2]) is to renormalize the tun- 
neling term. To restore the lattice symmetry broken 
by the time-dependent potential it is useful to intro- 
duce a vector potential A(i) = J^^ V{t')dt' and transform 
the solutions of the time-dependent Schrodinger equa- 
tion as = exp [-iA{t)x\'^{t). Note that the tunnel- 
ing terms, combination of creation and destruction op- 
erators on neighboring sites, acquire a non-trivial phase. 
On the other hand, the space dependent phase induced 
by the vector potential cancels out in the interaction and 
disorder terms, as they only depend on creation and de- 
struction operators on the same site. In the limit a; ^ oo 
one can average the phase over one period and show that 
the time-periodic external potential in Eq. ([3|) leads for 
to = to an effective time-independent Hamiltonian sim- 
ilar to Hq with an effective tunneling [2,, £, _19, 

Jeff = Je±^^/-^^"^Jo(^M, (9) 

where JoiV/uo) is the Bessel function of zero order. For 
finite but high-frequencies uo ^ U^J^ W ^ the Hamiltonian 
Hq with the approximate effective tunneling rate of Eq. 
([9]) still describes with high accuracy the cyclic states of 
the time periodic Hamiltonian ([2]) at to- For the zeroes 
of the Bessel function, the tunneling rate, although not 
exactly zero, is very much suppressed. Note that in this 
work we do exact time dependent calculations and do not 
use the effective Hamiltonian approximation. As shown 
below this approximation is very accurate and the main 
effect of the time-dependent potential Eq. (|3|) is to pro- 
duce a renormalization of the effective tunneling in the 
lattice. 

For an ordered lattice {W = 0) the Hamiltonian Eq. 
([2]) drives a transition [H, Q from an initial SF state into 
a MI state by adiabatically ramping the potential V to 
a value V/uo ^ 2.4, which corresponds to the first zero 
of the Bessel function, where Jeff = 0. In a disordered 
one-dimensional lattice with no interactions this oscilla- 
tory force can be used for increasing the localization of 
particles in the Anderson state in the regime of fast os- 
cillations [7], while for slow oscillations, it can be used to 
decrease the localization [22|, [23| . 

Here, we combine the oscillatory driving force, the ran- 
dom potential, and the interatomic interactions to ex- 
plore the possibility of localizing interacting particles not 
only in a MI state, but also in other insulating states. 
We explore in a controlled fashion the full parametric 
dependence of the ground state of interacting bosons in 
the presence of disorder. As we can effectively control 
the value of the effective tunneling matrix elements be- 
tween the lattice sites we can move between the different 
states in the diagonal lines of the corresponding phase di- 
agram U/J vs. W/J (i.e. we keep U /W constant while 
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FIG. 1: Quasienergy levels of the Floquet Hamiltonian as a 
function of the amplitude of the fast oscillatory force V for 
different values of the disorder potential strength a) VK = 
b) W/J = 0.1 and c) W/J = 1 for AT = L = 5, U/J = 0.5, 
uj/ J — 2b. For high frequencies the quasienergy levels e OL cor- 
respond to the eigenenergies of the static Hamiltonian in 
Eq.([T]) with effective tunneling Eq. (|9]). The structure of the 
quasienergy levels increases complexity for increasing values 
of disorder. The ground state remains gapped from the rest 
for all values of W allowing for the adiabatic following of the 
ground state. The properties of the ground state change dra- 
matically with V . Close to V /ijj — 2.4 the effective tunneling 
is suppressed and the corresponding ground state is localized. 

changing the effective U/J). For this control through ex- 
ternal periodic forces to be of any use, one should take 
into account the velocity of the ramping needed for the 
transition between states to take place. If one increases 
the amplitude of the external force adiabatically the sys- 
tem follows the eigenstates (cyclic states) of the Floquet 
Hamiltonian. We show that adiabatic behavior can be 
achieved for realistic parameters of our system. 

III. ADIABATIC DRIVING TO LOCALIZED 
STATES 

We start at t = from the ground state of Hq and 
switch on the time-dependent potential Eq. ([3]) assuming 
a discrete ramp Vm = AFm with m = 1, . . . , A/'g. For 
each Vm we solve the Floquet eigenvalue equation and 
obtain the corresponding quasienergies and cyclic states. 
Assuming that for each Vm we let the system evolve dur- 
ing a whole period T the evolved state reads 

I ^{t = N,T)) = U{T, 0; Vm)^ \ ^(0)) , (10) 

where U is given by Eq. ([8|). We have checked that in 
the adiabatic limit we obtain similar results for a contin- 
uous ramping of the external potential. According to the 
generalized adiabatic principle for Floquet Hamiltonians 
of finite size, we consider that the evolution is adiabatic 
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FIG. 2: (Color online) a) Quasienergy levels as a function of 
the amplitude of the fast oscillatory force V ior N = L = 5, 
U/J = 0.5, Lv/J = 25 and W/J = 0.1. At V/lv = 2.4, the 
effective tunneling is almost suppressed and the system local- 
izes. The quasienergy level structure shows a gapped struc- 
ture at this point, b) Evolution of the momentum distribu- 
tion of the state | ^(t)) when the fast oscillatory potential is 
ramped linearly with step AV/ J = 0.06. The evolution is adi- 
abatic as the state follows the ground energy eigenstate of the 
system. The momentum distribution changes from a peaked 
structure to a flat structure when the tunneling is suppressed. 

when the state | ^(t)) follows the same cyclic state during 
the ramping 

I {to, niAV) I U{{m - 1)T, 0) | ^{0)) | ~ 1 (11) 

for a particular a. As we will see below, the variation of 
the quasienergies with potential amplitude is such that, 
in general, the levels become closer when V/co ^ 2.4 
and the effective tunneling is coherently suppressed. 
In our calculations we have used only one ramping 
speed such that adiabatic evolution is assured in this 
region. Note, however, that the ramping could be faster 
initially for low values of V and slowed down at the end. 
Calculations have been performed using the full spatial 
basis of Ho and a truncated Fourier basis. The size of 
the Fourier basis is adjusted to assure convergence. We 
use the Arnoldi-Lanczos algorithm [23] to keep only the 
central eigenvectors with quasienergies Sa- 

To keep the problem computationally tractable we re- 
strict ourselves to systems of size L = 5. The results 
shown below are typical examples that demonstrate that 
the coherent control of our system is feasible. We dis- 
cuss the expected behavior for increasing system sizes 
and conclude that adiabatic following is expected to oc- 
cur also for realistic experimental sizes. Note that for 
small system sizes, some possible phases of the disor- 
dered BH model may not appear p^. Also, we want 
to be cautious when identifying phases which occur in 
many-particle systems [lo|, 13 ^^^^ ground states of 
the Hamiltonian Hq for small system sizes [12]. Our small 
system ground states are just precursor of the phases that 
appear in a large particle system. In this work we are in- 
terested in the time-dependent phenomena and our aim 
is to show that one can move adiabatically between the 
ground states of the effective BH Hamiltonian for varying 
values of the parameters U/J and W/ J with U/W fixed. 



This would correspond to an adiabatic transition along 
the diagonal lines of the phase diagram U/ J vs. W/ J for 
a many-particle system of larger size. 

All the results shown here correspond to parame- 
ters ^ = 0, and C = 0.002. To illustrate the 
effect of disorder in the system we present results 
for one realization of the random potential = 
{0.6804, -0.4182, -0.8932, 0.3686, -0.3248}. Any real- 
ization of disorder that completely breaks the symmetry 
of the system would yield similar adiabatic time scales 
although the precise value of the ramping time might de- 
pend on the exact form of the random potential. Specifi- 
cally the behavior of the time scale with occupation num- 
ber and disorder ratio W/U do not depend on the actual 
realization of the random potential. 

To characterize the state | ^(t)) we use the one-particle 
density matrix and the momentum distribution 

Pij = {^\alaj\^) (12) 
p, = J2e-"'^'-'H^\ala,\^) (13) 

(14) 

where Lk/{27r) = -M/2, ..,M/2. Both quantities evolve 
in time, as the fast oscillatory potential is ramped up. 

For an ordered system (W = 0) with commensurate 
filling, it has been shown theoretically [2] and experimen- 
tally [3] that the adiabatic ramping of the fast oscillatory 
potential in Eq.(j3|) can drive a SF-MI transition. In the 
presence of disorder, the reflection invar iance of the sys- 
tem is broken, the properties of the states are different 
and there is no guarantee that the adiabatic time scales 
are within experimental reach. We show in fig. [1] the 
increasing complexity of the quasienergy levels with dis- 
order strength. The breaking of symmetries makes the 
quasienergy spectrum more dense but strengthens the re- 
pulsion between energy levels. 



A. Weak disorder W/U < 1 

When the disorder is weak, in the limit U/J^O, the 
ground state of the system is predicted to be a MI [l2|, 
14], an insulating state with homogeneous distribution 
across the lattice and a gapped excitation spectrum. Our 
calculations show that one can drive a transition from 
an initial SF state at t = into a Ml-like state when 
we add small disorder W/U = 0.2 to the system. The 
quasienergy spectrum and the momentum distribution of 
the evolved state are shown in figure O The one-particle 
density matrix smoothly changes from a flat structure to 
a diagonal shape for V/co close to the first zero of the 
effective tunneling in Eq. (j9|) as shown in fig. [3] a) and 

b). 

Commensurate filling implies that the ground MI state 
at Jeff = is non-degenerate. As shown in fig. O the 
lowest quasienergy state is always gapped from the next 
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a) W/U=0.2 t=0 



b) W/U=0.2 t=t 




ramp 




a) W/J=0 J , =J 

' eff 




b) W/J=0 J^^^=0 




c) W/U=2 t=0 




d) W/U=2 t=t 



ramp 




c) W/J=0.1 J =J 

' eff 




d) W/J=0.1 J^^ =0 




e) W/U=140 t=0 




f) W/U=140 t=t 



ramp 




FIG. 4: (Color online) One-particle density matrix pij of the 
initial (left) state at t = and final states (right) at t = tramp 
with and without disorder for L = 5, = 4, cj/J = 25 and 
[//J = 0.5 . The initial states in figures a) and c) have very 
similar long-range coherence properties, while only the final 
state in d) with finite disorder W is localized. 



FIG. 3: (Color online) One-particle density matrix pij of the 
initial (left column) state when Jeff = J and final (right col- 
umn) state when Jeff = J for different values of the disorder 
strength and = L = 5. a) and b) cj/J = 25, WjJ = 0.1, 
AV/J = 0.06 c) and d) uj/J = 25, W/J - l,Ay/J = 0.01 e) 
and f) uj/J = 37, W/J = 7, AV/J = 0.3. 



excited states enabling the adiabatic transition. Remark- 
ably this behavior persists in the presence of weak disor- 
der and we obtain similar time scales for = and 
for W/U = 0.2 . Our calculations for TV = 5 yield 
AV/J = 0.06 which implies Ns = 2.4 x u/AV = 6000 
and total ramping times for V (tramp) = 2.4cc; given by 

AFtramp = 27r X 2.4. 

We thus obtain Jtramp ~ 250 which for typical J ~ kHz 
leads to experimentally achievable ramping times of the 
order of tenths of seconds. 

For non-commensurate fillings localization is not pos- 
sible for the proposed scheme for W = [IB]- We show 
in fig. mthe one-particle density matrix of the initial and 
final state for AT = 4 and W = and VI^//7 = 0.2. Figures 
m a) and c) show that both initial states have long-range 
coherence. The final state for = (fig. [Hb)) has long 
range coherence and indeed is a symmetrized combina- 
tion of all the Fock states with one particle per site and 
one hole. On the contrary, if we add disorder to the sys- 
tem, the ground state for Jeff = is the Fock state with 
one particle per site and the vacancy located in the site 
with the highest disorder as shown in fig. [4]d). Due to 
the presence of disorder, the degeneracy is thus lifted and 
one achieves adiabatic passage into this localized state for 
a ramping with potential step AV/J = 0.001. 




FIG. 5: (Color online) a) Quasienergy levels for N — L — 
5, U/J = 0.5, uj/J = 25, W/J ^ 1. b) Evolution of the 
momentum distribution of the ground state for a potential 
with a linear ramp of step AV/J — 0.01. 

Calculations with large particle numbers [TO, [lOl P^^" 
diet a BG phase between the SF and the MI states for 
intermediate values of J /U and W/U < 1. As shown in 
[l^ these BG states disappear when the system is small. 
We do not observe this intermediate state in our calcu- 
lations in agreement with the results in [12]. For larger 
systems where the intermediate BG phase appears, the 
fast oscillatory force Eq. (|3]) could be used to reach either 
MI states or BG states by just controlling the final value 
of the ramping potential. 



B. Intermediate disorder W/U ^ 1 

For stronger values of disorder, W/U ^ 1, the pre- 
dicted ground states of the disordered BH model for 
U / J ^ are BG states [H, [l^ • These states are insulat- 
ing, with a spatial distribution which follows the disorder 
potential pattern, and show a continuous spectrum. 
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Our time-dependent calculations show that we can 
drive an adiabatic transition from an initial SF-like state, 
the ground state of at t = 0, into a BG-like state at 
a final time when the effective tunneling Jeff ~ 0. The 
quasienergy spectrum and the momentum distribution of 
the evolved state are shown in figure [5l The one-particle 
density matrix smoothly changes from a flat structure 
peaked around the deepest well induced by disorder (see 
fig. [3]c)) into a diagonal shape for Vjijj close to the first 
zero of the Bessel function Eq. (j9|) as shown in fig. [3]d). 
To completely characterize the BG-like state obtained in 
this manner from the Ml-like state obtained in the pre- 
vious section for smaller values of the disorder we show 
in fig. [6] the spectrum obtained when the tunneling is 
coherently suppressed for Vjijj = 2.4. We observe the 
gapped MI spectrum for W jU = 0.2 and the continuous 
spectrum for the BG obtained with WjU = 2. 

We now look at the effect of the occupation number on 
the adiabatic time scale. For commensurate filling, the 
ground state which was always gapped from the excited 
states for small values of disorder, becomes now closer to 
the higher excited states and the time scale needed for 
adiabatic ramping increases. For the particular example 
we have calculated (A^ = 5 shown in fig. [5j), we need to 
go to larger time scales AV/J = 0.01 for W/U = 2 in 
order to observe adiabatic behavior. Note however in fig. 
[1] c) that the lowest energy state is still gapped across 
the transition. For non-commensurate filling, disorder 
breaks the degeneracy of all the Fock states with similar 
occupation patterns for the effective Hamiltonian with 
^eff = 0. Increasing disorder strength implies increasing 
gap and we thus obtain faster ramping scales. For N = 
4 we obtain AV/J = 0.01 for W/U = 2, one order of 
magnitude faster than for W/U = 0.2 and of the same 
order than in a system with N = 5 and W/U = 2. This is 
in contrast to the small disorder case in Sec lIII A| where 
the ramping depends strongly on the filling factor. 

We expect this behavior to persist even for larger sys- 
tem sizes because in a disordered potential the sites with 
minimum energy, ~ —W/2, are randomly distributed 
across the lattice. The first excited states correspond 
to particles distributed in those minimum energy sites 
which in principle are far away in the lattice and cannot 
be coupled to low orders by the nearest neighbor tunnel- 
ing term. 



C. Strong disorder W/U > 1 

In the regime where the disorder dominates, the con- 
densate state which appears for large values of J/U is 
not extended over the whole space, contrary to previous 
cases. For J/U large, the ground state is a condensate 
localized in the region of the lattice most distorted by 
the disorder. This localized condensed state which ap- 
pears in our small system resembles the properties of AG 
states appearing in larger systems [l2[, where particles 
condense into localized Anderson states with an expo- 




FIG. 6: (Color online) a) Quasienergy spectrum Sa (ot indi- 
cates the quasienergy ordering) of H(t) when V/uj = 2.4 such 
that Jeff = for W/J = 1 (A) and W/J = 0.1 (Q)- The 
continuous line for W/U = 2 resembles a BG-like behavior, 
while the discrete spectrum obtained for W/U = 0.2 resem- 
bles a Ml-like spectrum. {N = L = 5 and U/J = 0.5) b) 
Inset: Eigenvalues of the one-particle density matrix of the 
ground state of the initial Hamiltonian at t = (Hq) with 
U/J = 0.05,VK/J = 7 shown in fig. 3 e). Square of the most 
populated eigenfunction of the one-particle density matrix in 
fig. 3e). 




FIG. 7: (Color online) a) Quasienergy level structure for N = 
L = 5, U/J = 0.05, uj/J = 37, W/J = 7. b) Evolution of the 
momentum distribution of the ground state for a potential 
with a linear ramp of step AV/J — 0.3. 



nential decay distributed across the lattice in different 
islands. We calculate the one-particle density matrix of 
the initial ground state for W/U = 140 shown in fig. [3]e). 
The eigenvalues of the one-particle density matrix (inset 
of figlUb)) show that only one state is macroscopically 
occupied. This state, shown in figlUb), decays exponen- 
tially in space. The quasienergy spectrum and the mo- 
mentum distribution of the evolved state depicted in fig. 
[3 show that the system evolves from a state with finite 
condensate fraction into an insulating state. Note that 
for bigger system sizes, both the initial and final states 
are insulating, and our fast oscillatory force is able to 
drive a transition from states with off-diagonal terms in 
the one-particle density matrix around the deepest sites 
of the lattice into states with only diagonal terms in the 
one-particle density matrix, as shown in figiSl This fi- 
nal diagonal states appearing for Jeff ~ resemble BG 
insulating states. 

The one-particle density matrix shown in fig. [3] peaks 
around the deepest lattice site and looses its off-diagonal 
elements when the system becomes maximally localized. 
For = 5 we observe that the ground state remains al- 
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ways the ground cyclic state of the quasienergy structure 
shown in figl5]for ramps with /SV/ J = 0.3 which would 
mean tramp ^ 10~^s. In this regime when W/ J is large, 
the quasienergy spectrum in fig. [71 is dominated by disor- 
der for all values of the ramping potential V/ J. Even at 
t = when the tunneling is largest Jeff = J the hopping 
is mainly taking place around the deepest potential well. 

For larger system sizes, where the deepest sites are in 
principle far away, the ground state would be a prolifer- 
ation of islands of the type considered here. The near 
energy states appearing in the spectrum as the system 
size is increased are not connected by the nearest neigh- 
bor tunneling matrix term of the Hamiltonian and one 
can therefore expect that the AG-BG-type transition on 
each island takes place independently. We thus expect 
that the transition between ground states can be realized 
adiabatically within experimentally feasible time scales. 

D. System size considerations 

When the system size increases the quasienergy spec- 
trum becomes more and more dense. The requirement of 
high frequencies is that uo J^W and thus when N 

increases the quasienergies of the high energy states can 
cross that of the ground state; i.e. the frequency co may in 
principle induce resonances between different states. Yet 
the states which are coupled by uo have very different oc- 
cupation number patterns and do not have any effective 
overlap. This type of accidental crossings are traversed 
diabatically for any physically relevant time scale as dis- 
cussed in [2]. 

All the adiabatic time scales calculated in this work in- 
volve very small particle numbers. Due to computational 
restrictions, any exact time-dependent calculation for the 
BH model involves small number of particles. This limi- 
tation has not prevented small calculation predictions to 
be proven in real experiments with much larger number 
of particles [2|,l3l and no disorder. The disorder potential 
we consider is white noise of strength W which means 
that for a large system the average energy gap between 
neighboring sites (the ones coupled by the Hamiltonian) 
is of order W/2. In other words, for a large system, lattice 
sites distorted by —W/2 are very far on average. Thus, 
most energy levels which appear close to the ground state 
cannot be coupled to it to low orders indicating that our 



small system calculations can be extrapolated to larger 
system sizes. 

IV. CONCLUSIONS 

We have shown that the adiabatic ramping of a con- 
stant force rapidly oscillating in time allows for the com- 
plete scan of the parametric dependence of the ground 
state of a disordered BH model. As shown previously 
the inclusion of an ac field in the Bose-Hubbard Hamil- 
tonian is equivalent to a renormalization of the effective 
tunneling for the time independent spatial Hamiltonian. 
A succession of oscillatory fields where the amplitude is 
slowly increased gives rise in the adiabatic limit to a suc- 
cession of effective Hamiltonians from some initial J to 
the final Jeff = and thus allows for a complete scan of 
the parameters U/J and W/J with W/U constant. 

We have found that in spite of the increasing com- 
plexity of the quasienergy diagram with disorder, adi- 
abatic following of the ground state is feasible with a 
similar time scale as in a system without disorder. Adi- 
abatic following of the ground state enables the transi- 
tion from states with long-range coherence into localized 
states both for commensurate and non-commensurate 
fillings of the lattice potential. For strong disorder we 
have shown that the transition between insulating states 
localized in the most distorted regions of the lattice is 
possible. We have found that the addition of disorder 
shortens the time scale needed for adiabatic ramping if 
the disorder is weak or strong compared to the interac- 
tion strength. When the disorder is comparable to the 
interaction strength, adiabatic behavior is still feasible. 

Our results indicate that this type of transitions be- 
tween different regions of the phase diagram can be real- 
ized adiabatically for experimentally feasible time scales 
in systems with large number of particles. The adia- 
batic ramping of a fast oscillatory force could be used 
not only to control localization of disordered bosons but 
also to scan richer phase diagrams predicted for bosonic 
or fermionic systems in optical lattices. 
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